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Abstract 

We reexamine Liouville equation based turbulence theories proposed by Qian [Phys. Fluids 26 , 
2098 (1983)] and Edwards and McComb [J. Phys. A: Math. Gen. 2 , 157 (1969)], which are 
compatible with Kolmogorov spectrum. These theories obtained identical equation for spectral 
density q{k) and different results for damping coefficient. Qian proposed variational approach and 
Edwards and McComb proposed maximal entropy principle to obtain equation for the damping co¬ 
efficient. We show that assumptions used in these theories to obtain damping coefficient correspond 
to unphysical conditions. 
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I. INTRODUCTION 


Edwards |2j proposed turbulence theory based on Liouville equation for joint probability 
distribution function of Fourier modes ^^(k, f) of velocity field governed by forced Navier- 


10| were proposed 


Stokes equation. Following Edwards, a few more turbulence theories js, 4 ^ 
to solve Liouville equation and were reviewed by Leslie Q and McComb 8(|. Edwards theory 
seeks Fokker-Planck model representation for Liouville equation and obtains closed set of 
equations for spectral density q(k) and damping coefficient (total viscosity) cj(k) for sta¬ 
tionary, isotropic turbulence. The theory failed to be consistent with Kolmogorov spectrum 
1, 5| and the failure was attributed to equation for u(k) jsj. As a modification to Edwards’ 
theory, Edwards and McComb proposed principle of maximal entropy to derive an equa¬ 
tion for u(k) compatible with Kolmogorov spectrum. Within the Fokker-Planck framework 
of Edwards, Qian proposed variational approach and obtained different equation for 
damping coefficient consistent with Kolmogorov spectrum. Instead of using Ua(k,t), Qian 
used real dynamical modal variables X* and their governing equations which were suggested 
by Kraichnan and later utilized by Herring 4| for his self-consistent turbulence theory. 
Herring’s theory also uses Liouville equation and obtains equation for q(k) identical to equa¬ 
tion obtained by Edwards. In this paper, we reexamine Kolmogorov spectrum compatible 
theories proposed by Qian 101 and Edwards and McComb to obtain damping coefficient. 


We show in the next two sections that assumptions made in these theories correspond to 
unphysical conditions. 


II. VARIATIONAL APPROACH BY QIAN 

For discussion purpose, hereafter we refer to Qian’s variational approach as Q83. We 
use (Q83; to represent equation number (#) in Q83 paper (l0|. Qian based his theory 
on governing equation for real dynamical modal variables Xj for stationary, homogeneous, 
isotropic turbulence, written as 


dXi 

dt 


— (z/j — ^ AijmXjXmi (Q83; 8) 


( 1 ) 


where r'^X, represents external driving force. Einstein summation convention of repeated in¬ 
dices is not utilized in Eq. ([1]) and in this section. It should be noted that the real dynamical 
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modal variables and their equations were first suggested by Kraichnan within the context of 
hydromagnetic turbulence 6j. Qian’s theory seeks Langevin model representation for 
isotropic turbulence in the form 


—Xi ~ -niXi + fi, r]i = Ci + {i^- Vi) (Q83;12) 


( 2 ) 


by using 


AjjmXjXm — —QiXi + fi (Q83; 11) (3) 

j,m 

where —QXi is dynamical damping term and fi is white noise type forcing term. Qian pro¬ 
posed variational approach to obtain damping coefficient r]i. The approach yields equation 
for r]i by minimizing a function I{r]i), written as 


and using 



(Q83;29) 


— = 0. (Q83;28) 

dr]i 


(4) 

(5) 


Here ( ) represents ensemble average. Qian considered variation in / under constraint fi = 
constant, where fi is related to {Xf) by 


Xf) = 


Vi - u, 


Vi 


(Q83;25) 


( 6 ) 


Also, fi is proportional to spectral density q{k) [10 1. 

We now show that for I = Ifqf) and under the constraint of fi = constant, 



(7) 


within the framework of Langevin model considered by Qian. Consequently, the use of Eq. 
([5]) to obtain rji is in error. For stationary turbulence, solution of Eq. ([2]) suggests 


2,i(.y) = Fi 

in which correlation of white noise forcing term 


{mMt')) = F,d{t-t') 


( 8 ) 

(9) 
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is utilized. Here — is Dirac delta function. Using Eqs. ([3]) , (jl]) and ([9]), function /(r/i) 
can be written as 


^ = E ( E - i-QXt) ) = E (hh) = 5(0) ^ F.. 


( 10 ) 




Further, using Eqs 


), (ITOl) and for 0* = constant, we can write 


dl dF d2r]j (xf) 


This Eq. ffTTj) suggests that for all i 


dl 


drji 


T^O 


( 11 ) 


( 12 ) 


as 7^ 0. In view of this, use of Eq. ([5]), i.e. ^ = 0, in Q83 to obtain rji is in error and 
corresponds to unphysical condition = 0, Vi for stationary turbulence. 

Now we suggest possible modihcation within the framework of Q83. Consider a function 
V, written as 

_ _ 

(13) 




which satishes an exact condition 


j j 


dV 

dr]i 


(14) 


This condition along with Eqs. (0]) can be used, instead of Eq. ([5]), to obtain equation for 


Vi¬ 


lli. MAXIMAL ENTROPY PRINCIPLE BY EDWARDS AND MCCOMB 

For discussion purpose, hereafter we refer to Edwards and McComb’s theory as EM69. 
We use (M90; jf) to represent equation number [jf] in McComb’s book j^. Edwards and 
McComb Q considered stationary, homogeneous, isotropic, turbulence inside a cubic box of 
side L. Their Liouville equation based theory uses equation for Fourier modes ^^(k, f) of 
the velocity held Ua{ix.,t) governed by forced Navier-Stokes equation, written as 




u, 


,(k,f) = + /„(k,f), (M90;4.81) (15) 
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where fa represents external driving force, k is wavevector and k‘^ = |kp. The Einstein 
summation convention for repeated Greek indices is utilized while writing Eq. flT^ and 
in this section. Edwards and McComb js] theory seeks Fokker-Planck model equation for 
Liouville equation. The model equation contains two model parameters, namely r{k) and 
s{k), which account for contribution of nonlinear term in Eq. 0151) . The dynamical damping 
coefficient r{k) is related to damping coefficient u{k) by 


uj{k) = i/k‘^ + r{k). (M90;6.80) (16) 

The coefficient s{k) accounts for correlation of white noise forcing in the Langevin equa¬ 
tion for Fokker-Planck equation. Within the framework of Edwards and McComb and for 
stationary turbulence, ui{k) and s{k) are related by 

2u{k)q{k) = d{k) (M90;6.84) (17) 


where 

d{k) = W{k) + s{k) (M90;6.79) (18) 

and W{k) accounts for correlation of forcing term /^(k, t). The spectral density q{k) is 
defined by 

(t) = T>a/3(k)g(fc), (M90;6.85) (19) 

where Daisfk) = Sa 0 — EM69 considered entropy function S = S[q{k),u{k)] to obtain 
equation for oj{k) by maximizing S, corresponding to the condition 


Su{k) 


E 


ss 

Sq{j) 


SqU) 

5u{k) 


= 0. (M90;7.88) 


Here 


5q{j) d{k) 6 {k-j) 


+ 


6 d{j) 


(M90;7.89) 


( 20 ) 


( 21 ) 


Su{k) 2 u‘^{k) 2 ,u{j) Suj{k) 

and d{k — j) = 1 when k = j otherwise d{k — j) = 0. Edwards and McComb realized 
the difficulty in obtaining the second term on the right-hand side (rhs) of Eq. (|2T]) . After 
neglecting the second term, an approximate equation 

6 q{j) d{k) 6 {k-j) 


6 u;{k) 


2 u^{k) 


( 22 ) 
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was used for further calculation in EM69 


fl. 


This Eq. 


suggests that = 0 , Vj 7 ^ k. 


As a consequence, EM69 used following approximate equation 


5S 


d{k) 

2 u‘^{k) 


6 S 

Hk) 


= 0 (M90;7.90) 


(23) 


5uj{k) 

to obtain equation for (jj{k). 

We now show that the neglect of the second term on the rhs of Eq. fl2Tl) corresponds to 
unphysical condition. Consequently, the use of Eq. (|23l) to obtain u{k) is in error. Since 
EM69 is proposed for stationary turbulence. 


5 \ = Constant. 

j j 


\L) 


and from which we can write exact equation 


E 


5q{j) _ Hk) ^ ^ Sq{j) ^ ^ 


5u{k) Suj{k) 




6 uj{k) 


(24) 


(25) 


Substituting approximate Eq. 


of EM89 into Eq. 
d{k) q{k) 


and using Eq. 03. we obtain 

2 uj^{k) u{k) ^ ° 

and which is not correct for all k. In view of this, approximation used in EM89 to obtain 
u}{k) corresponds to unphysical condition q{k) = 0 and does not comply with conservation 
of energy Eq. (|2T1) for stationary turbulence where q{k) 7 ^ 0, Vfc. 

It should be noted that, within the framework of EM69, the unphysical behavior can be 
avoided if 


6 S 

duj{k) 


= 0 


(27) 


along with q{k) = constant, Wk is used instead of Eq. fl20|) . This means that S = S{u{k)) 
and second term on the rhs of Eq. (|20|) is equal to zero and is neglected. This kind of neglect 
by Qian in Q83 for function I{r]i) was considered mathematically incorrect by McComb j^. 
In our view, Eq. fl27l) can be considered as valid equation which seeks to optimize Entropy 
when energy of turbulence remains constant as q{k) = constant, V/c. 


IV. CONCLUDING REMARKS 


Within the Eulerian framework, a very few renormalized perturbation theories of turbu¬ 
lence are consistent with Kolmogorov spectrum 7|, |^. In this paper, we have reexamined 
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two such theories proposed by Qian 


10| and Edwards and Me Comb (3| and have revealed 


hidden unphysical conditions in these theories. We have suggested possible modifications, 
Eqs. (lT3j) . (11^ and fl27|) . to these theories but have not explored their usefulness in obtain¬ 
ing damping coefficient consistent with Kolmogorov spectrum. This will be explored in a 
broader context of our future work on turbulence theory development within the framework 
of Kraichnan’s direct interaction approximation j^. 
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